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We present the latest B hadron lifetimes, B® mixing, and D° mixing measurements using up to 3.0 fb -1 of 
data collected by CDF and DO experiments at Fermilab. The B® lifetime is measured from both J/"4><j> {CP 
admixture) and flavor specific channels, and the B c lifetime is obtained from semileptonic channels. Following 
the Bg oscillation frequency measurement at CDF in 2006, the DO collaboration now observes B® oscillation 
at a significance of about 3<r. Since the first D° mixing evidence established at B factories in 2007, CDF has 
observed D° mixing at a significance of about 4a level, the first time from hadron collider. 



1. Introduction 

Since the last FPCP conference in 2007, many in- 
teresting measurements have been carried at CDF 
and DO experiments thanks to the smoothly operating 
Tevatron. By March 30 2008, the Tevatron has deliv- 
ered about 4 fb _1 of integrated luminosity, while CDF 
and DO recorded about 3.2 and 3.4 fb _1 separately. 
The measurements related to these proceedings use 
data from 1.0 fir 1 up to 2.8 ftr 1 . The CDF and DO 
experiments are described in detail in Ref. 0, 0] • For 
these proceedings, the recent B® and B c lifetimes re- 
sults, B® and D° mixing parameters will be presented. 



2. Lifetimes measurement 

In the frame work of the Heavy Quark Expansion 
Theory (HQET), the inclusive decay rate of B hadrons 
is given by [3]: 
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in the limit of mj, — > oo or a free quark model, all the B 
hadrons should have same lifetime, given by the first 
term of Eq. [TJ The first correction arises from the 
kinetic and chromomagnetic operator which is at or- 
der of (Aqc£) /mb) 2 , the second correction arises from 
weak annihilation or Pauli interference which is at or- 
der of (A-QCD/mb) 3 ■ Both theoretical and experimen- 
tal uncertainties could be reduced if we measure the 
lifetime ratios of different B hadrons. Only the second 
correction will enter the lifetime ratio calculation. For 
example, Pauli interference will increase t(B + )/t(B°) 
and r(Af,)/r(i? ), and weak annihilation or scattering 
will reduce it. As a result, the lifetime ratios are not 
exactly unity for different B hadrons, and the precise 
measurements of the ratios can help test HQET at 
order of {Kqc d / m^) 3 . The current experimental sit- 
uation for tb+/tbo, tbq/tbo and T^/rga are shown 
in Fig. [TJ where one can see that the experimental re- 
sult and theoretical prediction agree well for the B + 



case. For the B® case, there is still some discrepancy 
between the prediction and measurements, while for 
the A;,, the theoretical prediction range is large and 
the experimental results have large uncertainty at this 
time. 
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Figure 1: Current lifetime ratios for different B hadrons 
with theoretical predictions from Ref. 0]. 



2.1. //' lifetime measurements 

The B° meson can decay directly or decay after it 
oscillates to B®. Time evolution of an arbitrary sys- 
tem a(t)\Bs) + b(t)\Bg) is governed by the Schrodinger 
Equation 



dt 



M -i- , , , 
2 ) \ b 



(2) 



where M and T are mass and decay matrices. Two 
mass eigenstates B\ and B° H appear as a result of 
this mixing property. The eigenstates have differ- 
ent decay widths and The average decay 
width is defined as Y = (r& + Yjj)/2, while the de- 
cay width difference is Ar = Yl — Yh- Ar also 
probes new physics since it's related to quantity <j) s 
by Ar = 2|r 12 |cos(0 s ), where S = arg(-M 12 /r 12 ), 
while Ti2 and M12 are matrix elements coming from 
the box diagram as shown in Fig. [2] New physics is 
likely to increase (j> s , so Ar could be smaller than the 
Standard Model prediction. 
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Figure 2: Feynman diagram for B® — B® mixing. 



2.1 .1 . Results from B° -> J/ip<p channel 

The data at CDF (1.7 ffcr 1 ) and at DO (2.8 ftr 1 ) 
were collected by di-muon triggers [J Q which pref- 
erentially selects events containing J/tp — » /J + /i~ de- 
cays. An artificial neural network is trained to sepa- 
rate signal and combinatorial backgrounds at CDF, 
which yields about 2500 signal events, while a cut 
based selection procedure is used to select data at DO, 
giving about 2000 signal events. The mass projections 
from CDF and DO are shown in Fig. [3] and Fig. SJ 
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Figure 3: The invariant mass distribution of B s candidates 
from CDF 

Since B° s is a pseudo-scalar meson, J/t/i and <fi are 
both vector mesons, the final states are admixtures of 
CP eigenstates, where S and D waves are CP even, 
P wave is CP odd. To separate the CP eigenstates, 
the angular distribution are used [5j. CP violation is 
predicted to be tiny in the Standard Model, so the 
CP phase is fixed to zero at CDF, while it's allowed 
to float at DO with flavor tagging. To extract lifetime 
t s and decay width difference Ar, the final fit is done 
with an un-binned maximum likelihood method on 
mass, lifetime and angular variables. 

The measured lifetime and decay width difference 
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Figure 4: The invariant mass distribution of B° candidates 
from DO 



results from CDF are @ 

t = 1.52 ± 0.04(stat) ± 0.02(syst) ps 
Ar = 0.08 ± 0.06(stat) ± O.Ol(syst) ps" 1 

where CP conservation is assumed, and the lifetime fit 
projection is shown in Fig. [5] At DO, the CP violation 
phase is allowed to float (with external constraint on 
strong phases from B° — > J/ipK*°), the results are 

r = 1.52 ± 0.05(stat) ± O.Ol(syst) ps 
Ar = 0.19 ± 0.07(stat)±g;g?(syst) ps" 1 

The lifetime fit projection is shown in Fig. [6] 
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Figure 5: Proper decay time projection of the B° — > J/ipcj) 
fit at CDF. CP even fraction is dominant, and the slope 
difference of even and odd curves indicates decay width 
difference. 

The systematic errors are controlled very well at 
both CDF and DO. Systematic errors considered at 
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Figure 6: Proper decay time projection of the B® — » J/ip(f) 
fit at DO. 
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Figure 7: K factor distribution for each channel. The DK 
channel is fully reconstructed, but treated in the same way, 
since the K track is misidentified and a K factor is needed. 



CDF include background angular distribution, mass 
model, lifetime resolution model, B° — > J/ifjK*° con- 
tamination, detector acceptance, and silicon detector 
alignment. Systematic errors considered at DO include 
procedure test, acceptance, reconstruction algorithm, 
background model, and detector alignment. 

2.1 .2. Results from B° -> D s ((I>ty)ttX channel 

This channel is different from B° B — > J/ip4>, be- 
cause it's flavor specific, i.e. can only decay into 
DJtt + X, while B® can only decay into D+ir~ X. By 
fitting the signal with one exponential function, the 
obtained lifetime is related to average decay width 
and decay width difference by 

r(S s °) /s = i(l + (^) 2 )/(l-(^) 2 ) (3) 

assuming no CP violation. So it can be used to con- 
strain r and Ar in the global fit. 

The result is obtained at CDF only, from both fully 
reconstructed B® — > D s ir and partially reconstructed 
channels such as B® — > D s p(tt + 7t°) where the tt° is not 
reconstructed. In the partially reconstructed chan- 
nels, a "K" factor is introduced to correct the proper 
decay time for missing mass and transverse momen- 
tum 



ct = 



L xy m 



rcc 
B 



K 



PT 



where L xy is the projection of decay length in the x — y 
plane, m r B c is reconstructed mass, pt is reconstructed 
transverse momentum, and 



K = 



m B PT 

m rec true 

m B p T 



where ms is the true mass and p$f ue is the true trans- 
verse momentum. The "K" factor distributions for 
different partially reconstructed channels arc obtained 
from Monte Carlo and shown in Fig. [7J 



The data (1.3 fb _1 ) were collected by displaced- 
vertex trigger 0] at CDF. The trigger takes the advan- 
tage of the long lifetime of B mesons and selects events 
within some impact parameter range (120 fim < do < 
1000 fim) with respect to the primary vertex. The 
trigger makes it possible to select B signal in large 
QCD backgrounds at CDF but also removes signal 
events which decay early. The lifetime bias is cor- 
rected using a trigger efficiency curve obtained from 
Monte Carlo. 

The lifetime is extracted from a two-step fit. First 
a mass-only fit is done to extract relative fractions of 
various modes and backgrounds. The corresponding 
mass shapes are obtained from either Monte Carlo or 
data. A lifetime-only fit is performed with fractions 
fixed from the previous mass-only fit, where the life- 
time is the only free parameter and the final result 
is [| 

t = 1.518 ± 0.041(stat) ± 0.025(syst) ps 

The lifetime results for control samples, such as B° — » 
D-(K+tt-tt-)tt+, B° -> D*-{D q (K+tt-)tt-)tt+, 
and B+ D (K+7r-)7r+ , all agree well with PDG 
values. 

Fig. [8] compares this result with all published re- 
sults from flavor specific channels. The previous PDG 
value is dominated by the DO result in 2006 from B® 
semileptonic channels. This new result agrees with 
the PDG 2007 value, but the central value is higher. 



2.2. B c lifetime measurement 

The doubly heavy B c meson is an interesting QCD 
laboratory, where both b and c quarks can decay as 
well as annihilate. The lifetime is expected to bc much 
shorter than light B mesons. Current theory predic- 
tions span from 0.47 ps to 0.59 ps [1(3], so experimental 
result with small uncertainty will bc useful to con- 
strain the theories. 
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The same lifetime measured at DO is hi 
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Figure 8: B® lifetime results from flavor specific channels. 



Both CDF and DO have measured the B c lifetime 
from B c — > J/iplX semileptonic channels, and data 
are collected by the di-muon trigger. At CDF both 
B c — > J/tpfiX and B c — > J/ipeX decays are recon- 
structed from 1.0 fb -1 data. At DO only B c -> J/ipfiX 
is reconstructed from 1.3 fb -1 data. Because of the 
missing momentum of the neutrino, a "K" factor is 
also needed from Monte Carlo to correct the recon- 
structed lifetime value. Here the mass is taken from 
the measured value [ll[ instead of the reconstructed 
value, so the "K" factor takes the form 



K = 



Pt(B c ) 



The main challenge of the analysis is the various 
and dominant backgrounds which include: 1) fake J/ip 
plus true lepton, 2) true J/ip plus fake lepton, 3) true 
J/ip and lepton from different b quarks, 4) prompt 
J/ip background, 5) residual conversion (for electron 
channel only). At CDF, the shape of the proper de- 
cay time distributions of all backgrounds are modeled 
and calibrated carefully from either data or Monte 
Carlo, and the lifetime is extracted from a lifetime- 
only fit. At DO, the mass shape of signal and back- 
grounds are also modeled, and the lifetime is extracted 
from a mass-lifetime simultaneous fit. 

The measured lifetime results at CDF from both 
muon and electron channels are [l2j 

ct = 179.1+27 2(stat) jLtm (muon channel) 

ct = 121.7j^!" jj(stat) jitm (electron channel) 

and the combined result from two channels is given by 

t(B c ) = 0.475±g;°,g(stat) ± 0.018(syst) ps 



t(B c ) = 0.448t°;o3g(stat) ± 0.032(syst) ps 

The measured B c lifetime values from CDF and DO 
agree with each other well, and since the Tevatron is 
currently the only place which can produce B c mesons, 
one can make a weighted world average with all the 
previous measurements from the Tevatron. The plot 
is shown in Fig. [51 and the weighted average is 



t{B c ) = 0.459 ± 0.037 ps 
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Figure 9: Comparison of B c lifetime with CDF and DO 
RunI results, and an weighted average (assuming no cor- 
relations among) those measurements. 



3. B° s Mixing 

As shown in Fig. ® oscillates through a box 
diagram. The oscillation frequency Am s is the mass 
difference of heavy and light eigenstates Am s = n%H — 
itll, which are related to the off diagonal element M\i 
of the mass matrix by 



Am s = 2\M 12 \ 



(4) 



With the oscillation frequency from the B° system, 
one can take the ratio to cancel most of the theoretical 
uncertainties and get 



Am, 



\V ts \ 2 



Am d \V td \ 2 



(5) 



which is very important for constraining the CKM 
unitary triangle with measurements of other CKM 
matrix elements and some theory inputs. 

The probability for a B® meson to mix or not as a 
function of time is proportional to 



P(t) B o_> B o b° oc 1 ± cos Am s t 



(G) 
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And one needs to identify the flavor of the meson at 
both production and decay time to see if it's mixed or 
unmixed. By reconstructing flavor specific channels, 
the flavor at decay time can be identical to the charge 
of the daughter particle. To identify the flavor at pro- 
duction time, two algorithms are generally used at the 
Tevatron. On the same side of the reconstructed B® 
meson, one looks at the charge of the kaon from frag- 
mentation processes. On the opposite side, one can 
look at the charge of leptons from semileptonic decays 
or the jet charge of the other b quark. However, many 
things can dilute the tagging power, for example, the 
kaon mis-identification and low efficiency on the same 
side, mixing of neutral B mesons and sequential b de- 
cays on the opposite side. Thus, a dilution factor T> 
is introduced, and it modifies Eq. [SJ 



P(t) 



B9^B?,B<> 



oc 1 ± T> cos Am, t 



(7) 



and (1 +T>)/2 gives correct tagging probability. 

Following the CDF measurement in 2006 with 
Am s = 17.77 ± O.lO(stat) ± 0.07(syst) ps, DO now 
has a new measurement with 2.4 fb _1 data. The first 
1.3 fb _1 belongs to Run Ha period and the other 
1.1 fb -1 belongs to Run lib when an additional sil- 
icon layer was installed at DO. Both hadronic and 
semileptonic decays are collected by inclusive single 
or di-muon triggers at DO. Most of the reconstructed 
channels and signal yields are listed in Tab. HI in ad- 
dition, 1.2 ffcr 1 of B° s -> fi+D-(JK°K-) decays from 
Run Ha give about 600 signal events. 



Table I Reconstructed channels and yields in B® mixing 
measurement 



Channel 


Run I la 


Run lib 


Total 




28238 ± 339 


16539 ± 239 


44777 ± 415 


e + Dj{^-) 


1142 ± 83 


548 ± 45 


1663 ± 102 


7r+D s -(07r-) 


159 ±13 


90 ±11 


249 ± 17 


li + D~{K*°K-) 


11649 ±661 


6449 ± 616 


18098 ± 903 



An amplitude A is introduced as shown in 

P{t)B°^B°,B oc 1 ± AD cos Am s t (8) 

where A can be fitted by fixing Am s at different 
points, the true Am s value is indicated when the fit- 
ted A is consistent with unity, otherwise a sensitivity 
at 95% C.L. can be defined by the probe Am s value 
at which 



1. 645(7,4 (Am se „ s ) = 1.0 



(9) 



Individual amplitude scans vs. probe Am s were done 
separately for different channels and combined using 
COMBOS program [1J|. The combined amplitude 
scan is shown in Fig. |TD] The likelihood profile is 



obtained from the combined amplitude scan using the 
formula UM 



ALnC = —LnC{Am s ) + LnC(oo) 



\ - A fit (Am s 



(10) 

and the profile is shown in Fig. [TTJ The measured 
Am s and its errors are derived by fitting a quadratic 
function to it, giving [l(| 



Am s = 18.53 ± 0.93(stat) ± 0.30(syst) 
which has a total significance of 2.9 a. 

D0 Runll Preliminary: [ L dt = 2.4 fb~' 

" " ii — m — r~r^ — i i i I i i i i I i J i i I i i i i I 



(11) 
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Figure 10: B® mixing combined amplitude scan. 



D0 Runll Preliminary: L dt - 2.4 fb 



stat error only 
total error 



A IrtL(min) s -4.0 
7 A mj, = 1 8.53 ± 0.93(statj : 0.30(syst) ps' 



12 14 16 18 20 22 24 26 28 30 

Am, (ps 1 ) 

Figure 11: B® combined likelihood profile, which is done 
with statistical error and total error on amplitude A sep- 
arately. An estimate of the Am, systematic error is ob- 
tained by subtracting, in quadrature, the error derived 
from statistical-error-only likelihood curve from that of the 
total-error curve. 



4. D° Mixing 

Since the discovery of charm quark in 1974, physi- 
cists have been trying to observe oscillations of the 
neutral charm meson. Unlike the kaon or B sys- 
tems, the oscillation frequency in the charm system 
is predicted to be very small. In the Standard Model, 
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D° mixing occurs through two processes. The "long 
range" mixing can be estimated using strong interac- 
tion models. The "short range" process is negligible 
in the standard model, however, exotic weakly inter- 
acting particle from new physics could enhance the 
short-range process. 

Recent D° mixing evidence has been found at the 
B factories in 2007. The Belle Collaboration found 
direct evidence by comparing the decay time distri- 
butions for D° decays to the CP-eigenstates K + K~ 
and 7r + 7r~ with that for the CP-mixed state K~tt + . 
The evidence found in the BABAR experiment is the 
decay rate difference of doubly-Cabibbo-suppressed 
(DCS) D° -> K+n- decay and Cabibbo favored (CF) 
D° — * K~ir + decay. The ratio of decay rates can be 
approximated as a simple quadratic function of proper 
decay time and mean D° lifetime, with the assumption 
of CP conservation and a small value of x — Am/T 
and y = AT/2T, 



R{t/r) = R D 



x 12 + y' 2 



(t/rf (12) 



where Rjj is the square of the ratio of DCS to CF 
amplitudes, x' and y' are linear combinations of x and 
V as 



V 



x cos 5 + y sin 5 
-x sin 6 + ycosS 



where 8 is the strong phase difference between the 
DCS and CF amplitudes. In the absence of mixing, 
x 1 and y' are both zero. 

At CDF, the same kind of evidence is found as in 
BABAR, but probed over a much wider D° decay 
time range. Events are selected with the displaced- 
vertex trigger using about 1.5 fb _1 data. The "right- 
sign" (RS) CF decay chain D*+ -> tt+P , P° -> 
K~tt + , and the "wrong sign" (WS) decay chain 



D 



*+ 



P u , P u 



K + tt are reconstructed. The 



relative charges of the two pions gives the "sign" of 
the decay chain and extra cuts are applied to reduce 
background to the WS signal from RS decays. The 
data are divided into 20 bins of t/r ranging from 0.75 
to 10.0 with bin size increasing from 0.25 to 2.0 to 
reduce statistical uncertainty at larger times. The ra- 
tio R is determined at each bin and a least-squares 
parabolic fit of the data is done with Eq. [T2J The 
fit is shown in Fig. 1121 and final results are shown in 

Table, nn In |- 
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Figure 12: Ratio of prompt D* "wrong sign" to "right 
sign" decays as a function of normalized proper decay time. 
The dashed curve is from a least-squares parabolic fit, the 
dotted line is the fit assuming no mixing 




Figure 13: Bayesian probability contours in the x 12 — y' pa- 
rameter space corresponding to one through four equiva- 
lent Gaussian standard deviations. The closed circle shows 
the unconstrained fit, the open diamond shows physically 
allowed fit (x' 2 > 0), the cross shows the no- mixing point. 



To determine the consistency of the data with the 
no-mixing hypothesis, Bayesian contours are com- 
puted with likelihood and a flat prior. The contours 
are insensitive to modest changes in the prior, and 
shown in Fig. 1131 The no-mixing point lies on the 
contour which excludes the region containing a prob- 
ability of 1.5 x 10~ 4 , equivalent to 3.8 Gaussian stan- 
dard deviations. Alternative procedures are also used 
to determine the probability for no mixing, and all 
give consistent results. 



Table II Fit results for the R(t/r) distribution 

Fit type I R D (1Q- 3 ) | y'(10~ 3 ) | a' 2 (10~ 3 ) \ x 2 /d.o.J 

Unconstrained 3.04 ± 0.55 8.5 ± 7.6 -0.12 ± 0.35 19.2/17 

Physically allowed 3.22 ± 0.23 6.0 ± 1.4 19.3/18 

No mixing 4.15 ± 0.10 36.8/19 



5. Summary 

The most precise B® lifetime and decay width dif- 
ference have been measured directly from the B® — > 
J/4"P channel from both CDF and DO experiments. 
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The measured lifetime is also consistent with the value 
from B® flavor specific channel. Together, the recent 
results confirm the theory prediction that Ts/t^ ~ 1. 
The lifetime measurement result for the heavy B c me- 
son can provide good experimental input for the the- 
ory calculation which still has a huge uncertainty. The 
fast oscillation of B° s mesons has been measured at 
CDF in 2006 and is now verified by the DO experi- 
ment with results from combined channels. While D° 
mixing is predicted to be quite small in the Standard 
Model, the evidence has been found at CDF followed 
the first observation at B factories last year. 
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